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ABSTRACT 

The importance of Partial Redistribution (PRD) in the modelling of the Lyman a 
and Lyman j3 emission lines of hydrogen in stellar atmospheres is examined using 
simple atmospheric models of a range of late-type stars. These models represent the 
sub-giant Procyon (F5 IV- V), and the two giants (3 Gem (KO III) and a Tau (K5 III). 
These stars are selected to span a wide range of surface gravities: 1.25 < logq < 4.00. 
The calculations are perfor med using the computer code MULTI ( Carlsson 1986| ) with 



the modifications made by Hubeny fc Lites (1995). It is found that PRD effects are 



highly significant, both in the direct prediction of the Lyman line profiles and in the 
application of hydrostatic equilibrium (HSE) to calculate the atmospheric electron 
density in static atmospheric models. 
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1 INTRODUCTION 

The modelling of the Lyman lines of hydrogen is extremely 
important for the study of stellar atmospheres. The Ly- 
1 man a line is the strongest single line in the ultraviolet 
spectra of cool stars and represents a significant energy loss 
mechanism from the high chromosphere and low transition 
region in most stars. The strength and width of the Ly- 
man lines, particularly Lyman a and Lyman f3, is such that 
they are import ant sources for fluore scence in the lines of 
other elements (Jordan & Judge 1984). This fluorescence is 



elled using the more complete theory of partial (frequency) 
redistribution (PRD). 



of particular importance in low gravit y stars such as a Tau 
(McMurry, Jordan & Carpenter |l999| ) . 

Unfortunately, the Lyman lines are difficult to model 
since the extreme optical thickness of the lines and the low 
electron densities in the region of line formation makes the 
treatment of photon scattering very important. The assump- 
tion of complete (frequency) redistribution (CRD), which 
is often used in the modelling of the lines of other atomic 
species, is unacceptable because most scattering events in 
a transition such as Lyman a are coherent (in a giant star 
such as f3 Gem in a region where the electron temperature 
T e ~ 10 4 K and the electron density n e ~ 10 9 cm -3 , only 
~ 1CP 4 of the scattering events are incoherent). This co- 
herence leads to a correlation between the absorbed and 
emitted frequencies in a scattering event and so introduces 
an explicit coupling of the radiation field at different fre- 
quencies. This significantly complicates the problem that is 
to be solved. Thus scattering in these lines must be mod- 
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Linsky (1985) gives a review of the work on PRD in 
stellar atmospheres. Early work on the Lyman a line in the 
solar chromosphere gave cl ear indications that PRD effects 
would be very import ant (Milkcy fc Mihalas 1973 ). How- 
ever, the treatment of Milkey fc Mihalas (1973)| is not en- 
tirely satisfactory because it is based on th e redi stribution 
function theory of Omont, Smith & Cooper (1972). This re- 
distribution function is widely used in the study of PRD in 
lines of Ca II and Mg II, but is not suitable for hydrogen 
because it assumes that the coherence fraction is frequency 
independent. However, the incoherence fraction is only fre- 
quency independent in the so-called Impact Limit of colli- 
sional broadening which is valid only in the very core of 
Lyman a. An adequate treatment of the line wings requires 
a more sophisticated atomic theory (under the conditions 
that exist in a giant star such as f3 Gem the Impact Limit 
is only valid for wavelength differences from line centre of 
< 10 -4 A). Such a theory exists in the form of the U nified 
Theory of Collisions (Smith, Cooper & Vidal lf969|). This 
theory is developed by Vidal, Cooper & Smith (1970) and 
specifically used to s tudy hydrogen in stellar atmosph eres b y 
Yclnik et al. (1981)| and Cooper, Ballagh & Hubeny (^98§). 
An alternative theory of collisi onal b roadening is discussed 
by Smith, Cooper & Roszman (1973). 



The application of the Unified Theory in the study of 
PRD in hydrogen has been carried out by several authors. 
Basri et al. (1979) used a Unified Theory treatment of both 
ion and electron collisions to compute a Lyman a profile 
to compare with solar observations. Their calculations were 
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based on the solar model of Vernazza, Avrett & Loeser 



(1973) which they modified to obtain close agreement with 



the observed profile. In their extensive study of solar chro- 
mospheric modelling, Vernazza, Avrett & Loeser (1981) also 
performed PRD calculations for hydrogen. In that paper, the 
authors experimented by setting different maximum coher- 
ence fractions and they found that enforcing a maximum 
coherence fraction of 98 per cent gave the best agreement 
with the observed data. However, they gave no physical rea- 
son why the incoheren ce fraction should be larger than the 



theoretical prediction. Luttermoser fc Johnson (1992) 



a PRD treatment of the Lyman a line when considering 
their joint photospheric/chromospheric models of the stars 
g Her and TX Psc. They find that a PRD treatment of 
hydrogen (compared to a CRD treatment) makes a signifi- 
cant difference to the proton density in certain parts of their 
atmospheric models and is t herefore significant to the ion- 
ization balance in hydrogen. (Jitenbroek (1989) introduced 
PRD to the radiative tran sfer code MULTI , whi ch in its orig- 
inal form is discussed by |Carlsson (1986) and Scharmer & 



Carlsso l (1985a, b) , The treatment of PRD used by (Jiten- 



broek (1989) was based on Impact Theory and is suitable 
for the Ca H H and K lines but not for hydrogen. More re- 
cently, [Hubeny fc Lites (1995)| modified the MULTI code 



extensively so that it could include the frequency dependent 
coherence fraction required for hydrogen. They discussed 
the modifications to the code and the theory behind them 
and also successfully applied their code to the solar chromo- 



spheric model C of Vernazza et al. (1981) 



The purpose of this paper is to present the results 
of appl ying the version of MULTI written by Hubeny & 
Lites (1995) to models of stars other than the Sun, so that 



PRD effects can be investigated in lower surface gravity and 
their importance to the modelling techniques can be studied. 
Several limitati ons remain in the approach used by Hubeny 



& Lites (1995) and these remain in the work here. The pri- 
mary limitation of the computer code is that it implements 
an angle averaged redistribution function. This means that 
it can only be applied realistically to a static atmospheric 
model and so all the models discussed here are static and 
contain no net flows. There is also a minor assumption, 
made for computational ease, that when calculating the col- 
lisional broadening, it is assumed that the proton density 
is equal to the electron density. In a typical high chromo- 
sphere/transition region this assumption is reasonable and 
so should not be a source of significant error. 



2 THEORY 

The detailed theory of PRD in hydrogen is discussed by 
Cooper et al. (1989)1 The theory they des cribe, which is 
that implemented by pHubeny fc Lites (1995) , includes an ac- 



ceptable treatment of resonant scattering in both Lyman a 
and Lyman f3 and also includes cross-redistribution from 
Balmer a to Lyman f3. The deta ils of the method are dis- 



cussed by [Hubeny fc Lites (1995) . A brief summary of those 



parts of the theory which are of particular relevance are 
given below. 

The study of PRD requires the calculation of the re- 
distribution function R which gives the probability that a 
photon of frequency v' will be absorbed and re-emitted as a 



photon of frequency v. It is convenient to work in the rest 
frame of the absorbing atom (on the understanding that, at 
a later stage an average over a Maxwellian velocity distri- 
bution of such atoms will be performed). Thus the required 
quantity is 

f{v')p{v',v) 

where / is the absorption profile (which describes the proba- 
bility that a photon of frequency v' is absorbed) and p is the 
probability that, given an absorption at frequency v' there 
is a subsequent re-emission at frequency v. 

The quantities / and p can be calculated from the quan- 
tum mechanical theory of radiation. As an example, for Ly- 
man a they are given by 



f ls 2p(Auj) 



and 



r 2p (Aoj) 



7T Aur + V 2p {Auj) 



(1) 



p ls2p (Au>', Alv) = (1- Ai s2p (Aa/)) 5(Aui- Aw') 

+Ai s2 p(Aw')/i s 2p(Aw) (2) 

where Auj and Aui' are the angular frequencies (measured 
relative to line centre) of the emitted and absorbed photons, 
T 2 p is the total half width of the 2p state and Ai S 2p is the 
incoherence fraction. 5 is the Dirac-delta function. 
The total half width is given by 



r 2p (Aw) = 72 P (Aw) + r 2p i s /2 



(3) 



where r2pi s is the radiative decay rate in Lyman a and 
72 P (Aoj) is the half width due to collisional broadening (for 
clarity, other broadening mechanisms are neglected here) 



Expressions for 72p(Aw) are given by Basri et al. (1979) 
The incoherence fraction is given by 



Al s 2p(A(j) 



72 P (Ao;) 
r 2p (Ao;) 



(4) 



The incoherence fraction is the probability that an atom 
will undergo an elastic collision (which randomly changes 
the energy and thus destroys coherence) before it decays. 
Thus Equation (2) means that a fraction A of scattering 
processes will undergo complete redistribution and a frac- 
tion 1 — A will undergo no redistribution. The incoherence 
fraction is clearly the physical quantity which determines 
the importance of PRD effects. It always lies between and 
1. If is it then the transition in question exhibits pure 
Rayleigh (coherent) scattering, while if it is 1 the system 
displays complete redistribution (CRD). Thus PRD effects 
are only likely to be important for values of the incoherence 
fraction that are <C 1. 

The incoherence fraction is determined by the broaden- 
ing paxam^tej^^^^Atj) . Th is quantity is discusse d in detail 
by [Yelnik et al. (1981) and |Cooper et al. (1989)|. It is obvi- 



ous that 72 P (Aoj) should be density dependent, since a lower 
density will lead to fewer collisions and hence a smaller in- 
coherence fraction. Therefore, it is to be expected that PRD 
should be significantly more important in a giant star (such 
as P Gem which has a surface gravity of 5.6 x 10 2 cm s~ 2 ) 
than in a main sequence star like the Sun (surface gravity 
2.75 x 10 4 cm s -2 ). It is for this reason that it is of particu- 
lar interest to examine the trend of PRD effects in different 
types of star. 
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The treatment of the Lyman (3 line is further compli- 
cated by the possibility of resonant Raman scattering involv- 
ing the Balmer a line. When a Balmer a photon is absorbed, 
exciting the n = 3 level, the n — 3 level may subsequently 
decay to n — 1 emitting a Lyman /3 photon. The frequency 
of the emitted Lyman j3 photon will be coupled to the fre- 
quency of the absorbed photon, unless a collision occurs to 
redistribute the n = 3 electron energy. A reliable treatment 
of the Lyman /? line must therefore account not only for 
significant coherence in resonant scattering in the Lyman j3 
line, but also for the effects of this so-called cross redistri- 
bution between Balmer a and Lyman j3. Equations (47) to 
(49) in Hubeny & Lites (1995) show how cross redistribution 



appears in the emission coefficient of Lyman j3. 

An additional problem occurs in the application of the 
Unified Theory to hydrogen - the Z-degeneracy of the en- 
ergy levels means that the Isolated Line approxi mation is 
not val id. A more general theory is de velope d by 



Burnett 



et al. ( p.980)| and Burnett & Cooper (|l980aj,b). This pre- 
dicts corrections to the Unified theory which take the form 
of additional correlation terms in the emission coefficients. 
These terms represent the collisional mixing of /-degenerate 
levels during far wing emission and are discussed in t he con- 

These 



(1989) 



text of stellar atmospheres by Cooper et al. 
terms are known to cancel out if the Z-degenerate levels are 
populated according to their statistical weights. It is known 
that in solar models the Z-degenerate levels are populated 
according to their statistical weights in the chromosphere 
and so these co rrelations terms were n ot included in the 
study made by Hubeny & Lites (1995). These correlation 



terms have also been discarded in the present work, but it is 
noted that, in principle, they may become more significant 
in the lower gravity stars considered here since the collision 
rates between the Z-degenerate levels will be lower due to 
the lower chromospheric particle densities. Estimates of the 
departure from population of Z-degenerate levels according 
to statistical weights have been made and are discussed in 
Section 4. 



3 MODELS 

PRD calculations have been performed using atmospheric 
models of three different stars, the low gravity giant a Tau 
(K5 III), the higher gravity giant /3 Gem (KO III) and the 
sub-giant Procyon (F5 IV- V). Basic parameters for these 
stars are given in Table pi. 



The a Tau model adopted is t hat of |McMurry (1999) 



This model was developed from the Kelch et al. (1978) chro- 
mospheric model and an observed emission measure distri- 
bution for the transition region. It was contructed using an 
approximate PRD approach, based on an truncation of the 
Lyman a profile. The properties of the high temperature 
material in a Tau are not well constrained observationally. 
The atmospheric model only extends out to a column mass 
density of 10 _5 g cm -2 and a temperature Te = 10 5 K which 
means that the optical depth in both Lyman a and Ly- 
man j3 is still significant at the top point in the atmosphere. 
When using this model to predict observations, the atmo- 
sphere is extrapolated outwards to zero optical depth, as- 
suming a constant source function. This estimate of the be- 
haviour of the atmosphere above the top of the model lim- 



Table 1. Fundamental stellar parameters. 



Parameter 


a Tau 


Gem 


Procyon 


Catalogue No. 


HD 29139 


HD 62509 


HD 61421 


Spectral Type 


K5 III 


KO III 


F5 IV-V 




3920 K a 


4865 K e 


6500 K9 


log g* 


1.25 6 


2.75 e 


4.00 h 


Radius 


44.3 R Q C 


8.93 R e f 


2.1 H ! 


Dist. from Sun 


19.96 pc d 


10.3 pc d 


3.5 pc l 



Blackwcll, Lvnas-Gray & Petford ( 1991 ) 



Bonncll fc Bell (1993) 



c Calculated using Blackwcll et al. (1991) 
d From Hipparcos CESA 1997). 



Drake fc Smith (1991' 



E 



Calculated using Mozurkewich et al. (1991) 
B kyrcs et al. (T974)( 



rake fc Laming ( 1995) 



Wood et al. (1996) 



its the reliability of the calculated line core profiles, up to 
the wavelength difference from line centre where the optical 
depth at the top point in the atmosphere becomes negligi- 
ble, which occurs at a detuning of around 0.4A for Lyman 
a and 0.3A for Lyman /3. This extrapolation is significant 
in the calculated integrated line fluxes, particularly the Ly- 
man /3 CRD flux which is dominated by the flux in the 
optically thick line core. The a Tau model does not account 
for the stellar wind which will contain a significant fraction 
of neutral hydrogen. The wind will introduce an additional 
source of optical depth above the top of the atmospheric 
model which further limits the physical reality of the upper 
boundary conditions in the model. Recent Far Ultraviolet 
Spectroscopic Explorer (FUSE) observations of a Tau will 
help to constrain the properties of the upper atmosphere of 
the star. 

The P Gem and Procyon models are new and have been 
constructed by the author in a similar fashion to the a Tau 
mode l. The models are based on the chromospheric mod- 
els of |Kelch et al. (1978)| and Ayres, Linsky & Shine (]l974|) 



respectively, and observed emission measure distributions. 
These new models are initially constructed using the as- 
sumption of CRD. The details of the construction of these 
models will be discussed in future papers. The new models 
extend to coronae with logT c = 6.3. The Lyman lines are 
optically thin at the top of these coronal models and so the 
accuracy of the calculated line cores are not limited by an 
extrapolation in the same way as in a Tau. The electron 
temperature as a function of the column mass density is 
shown for each of the models in Figure |l[ 



4 RESULTS AND DISCUSSION 

The Hubeny & Lites version of the MULTI code has been 
applied to each of the three atmospheric models described in 
Section |l| Each model has been run for four different cases, 



the same cases that were studied by Hubeny & Lites (1995) 
namely: 

(a) CRD 
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-6 



-5 



7F 



a Tau 



log column mass density 



(S Gem 



-4-2 2 

log column mass density 



Procyon : 



log column mass density 



Figure 1. The atmospheric models. In each case the electron 
temperature is given in K and the column mass density is given 
in g cm -2 . 



(b) PRD resonance scattering in Lyman a, 

(c) PRD resonance scattering in Lyman a and Lyman f3, 

(d) PRD resonance scattering in Lyman a and Lyman /3, 
and cross-redistribution between Balmer a and Lyman /3. 

For all twelve of these calculations, the physical depth 
and frequency dependent coherence fractions were used. 
Balmer a is treated in CRD throughout (Hubeny & 



Hcinzel 1984 show that PRD is not significant in determin- 



ing the line profile for a subordinate line, provided that the 
damping parameters for the upper and lower levels are not 
too different). The calculations were performed using a nine- 
level atomic model of hydrogen (t he same hydrogen model 
which was used by McMurry 1999). This model includes 28 
bound-bound transitions and 8 free-bound transitions. The 
Lyman continuum is calculated in detail while the others 
are treated using radiation temperatures. There is no fine 
structure in the atomic model (the fine structure splitting 
between components is much smaller than the width of the 
Lyman lines and is not observable). 

As discussed in Section 2, the PRD treatment employed 
here neglects certain correlation terms, the importance of 
which depends upon how closely the Z-degenerate levels are 
populated in accordance with their statistical weights (if the 
Z-degenerate levels are populated exactly according to their 
statistical weights these correlation terms have no effect). 
In order to attempt to estimate the relative populations of 
the /-degenerate levels basic statistical equilibrium calcula- 
tions have been applied to a simplified fine structure model 
of hydrogen. This model consists of only six levels (Is, 2s, 
2p, 3s, 3p and 3d). The transition rates between these levels 
due to spontaneous emission, stimulated emission, radiative 
absorption and electron collisions (both elastic and inelas- 
tic) were calculated for an electron temperature T e = 10 4 K 
(and the appropriate electron density for each of the three 
stellar models) using the radiation field from MULTI calcu- 
lations, the estimates for the elastic collision rates given by 
Cooper e t L jiL_(1 989)| and the inelas tic collision rates calcu- 



lated by 
Avrett 



Giovanardi fc Palla (1989) (note that while Chang 



Loeser 



1991 



show that there are inconsistencies in 



the rates calculated by Giovanardi & Palla 1989 involving 



states with n > 4 they do not find significant problems with 
the fine structure rates between states with n < 3). Under 
the conditions used for each of the three stars it was found 
that spontaneous emission, radiative absorption and elastic 
collisions are the important processes while stimulated emis- 
sion and inelastic collisions are not important under these 
conditions. The calculated relative level populations for each 
of the fine structure levels in n = 2 and n = 3 are given in 
Table ^. The ratios of statistical weights for the fine struc- 
ture levels are ui2 P /uJ2s = 3, (jJ3 P /oJ3 S = 3 and W3d/w3 S = 5. 
The results shown in Table ^ suggest that the fine structure 
level population ratios are likely to be very close to the ra- 
tios of the statistical weights (the greatest variation being 
4 per cent in the ratio of the 3d to 3s populations) and so 
the correlation terms may be safely neglected under condi- 
tions typical of the region where PRD effects are important 
in these stars. When considering these results, it must be 
borne in mind that the fine structure model that has been 
used is very simple: transitions between levels with principle 
quantum number greater than n = 3 have been ignored. It 
is estimated that, for example, radiative decays from n — 4 
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Table 2. Relative level populations for fine structure levels. The 
notation n n ; denotes the population of fine structure level nl. 



Star 


«.2p/«2s 


™3p/™3s 


™3d/™3s 


a Tau 


2.97 


3.09 


4.95 


P Gem 


2.99 


3.03 


4.83 


Procyon 


3.00 


3.00 


4.94 



to n — 3 will populate the n = 3 at a rate of around ten per 
cent of the rate at which n = 3 is populated by radiative 
absorption from n — 2 (which is the dominant excitation 
mechanism for n — 3 under the conditions in the simple fine 
structure model). 




1 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

AX 



4.1 Line Profiles 

The calculated Lyman a profiles for the four cases of each 
of the models are shown in Figure |^ and the corresponding 
Lyman j3 profiles are shown in Figure |§| These plots confirm 
that the treatment of redistribution has a significant effect 
upon the calculated hydrogen line profiles. Most of the dif- 
ferences between PRD and CRD results is due to the effect 
of coherent scattering (in the atom frame) in Lyman a which 
is included by PRD but not by CRD. It is clear that the pri- 
mary effect of a PRD treatment rather than a CRD treat- 
ment is a significant reduction in the strength of the Lyman 
a line wings. This effect occurs because the low incoherence 
fraction in PRD greatly reduces the number of photons that 
are scattered into the line wings compared to complete in- 
coherence in a CRD treatment. In all three stars considered, 
the treatment of Lyman (3 does not have a significant effect 
on the Lyman a profile (for each star cases (b), (c) and (d) 
give identical Lyman a profiles), however the treatment of 
Lyman a does strongly effect the calculated Lyman (3 profile. 
Coherent scattering (in the atom frame) in Lyman a affects 
the calculated Lyman ft as strongly as coherent scattering 
in Lyman (3 - coherent scattering in Lyman a traps photons 
in the Lyman a line which pumps the n = 2 level popu- 
lation above its CRD value, this increases the n = 3 level 
population since radiative excitation from n = 2 to n — 3 
(absorption of Balmer a photons) is a significant population 
process for n — 3. This increase in the n = 3 population 
leads to an increased flux in the Lyman f3 line wings. This 
effect is evident in all the case (a) Lyman j3 profiles. The 
effect of including PRD effects in the treatment of resonant 
scattering in Lyman /3 is to reduce the near wing flux by 
a similar mechanism to that which operates in Lyman a 
(see above). Finally, the inclusion of cross-redistribution in- 
creases the line w ing flux again. These eff ects were all noted 
and discussed by Hubeny & Lites (1995) for the solar case. 

The most striking result from Figures ^ and ^ is the 
significance of the trend with surface gravity. It is clear that 
the importance of PRD effects grows rapidly as stars of lower 
surface gravity (and hence density) are examined. This is 
to be expected because the elastic collision rates (and hence 
the incoherence fraction) are smaller in low density stars. As 
noted in Section |^, PRD is most different from CRD when 
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Figure 2. Calculated Lyman a profiles. The dotted line is case 
(a), the broken line is case (b), the dot-dash line is case (c) and 
the solid line is case (d) (some cases are indistinguishable: see 
text). All fluxes are in ergs cm _2 s _1 A — 1 at the stellar surface. 
AA is the wavelength relative to line centre in A. 
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Figure 3. Calculated Lyman j3 profiles. The dotted line is case 
(a), the broken line is case (b), the dot-dash line is case (c) and 
the solid line is case (d). All fluxes are in ergs cm~ 2 s _1 A _1 at 
the stellar surface. AA is the wavelength relative to line centre in 



the incoherence fraction is small, so it is to be expected 
that the greatest difference occurs in low surface gravity 
stars. (5 Gem and, in particular, a Tau show very much 
stronger PR D effects that Procyon or the Sun (Hubeny & 
Lites 1995). In both cases the reduction in the width of 



the Lyman a line is extremely significant and in a Tau the 
enhancement of the near wings of Lyman (3 is sufficient to 
make them stronger than the line core. Notice that even 
in these extreme examples of the importance of PRD, its 
inclusion in the treatment of Lyman (3 is not as significant 
as in the treatment of Lyman a. However, the effects on 
Lyman (3 are still significant and should not be disregarded. 

For comparison the models have also been run in 
MULT1 using two approximate PRD techniques. The first 
approximate technique (which will be referred to as case (e)) 
is the commonly used method of simulating PRD effects 
by performing CRD calculations using a truncated Voigt 
absorption profile. In the calculations presented here the 
absorption profiles for both Lyman a and Lyman /3 have 
been truncated at A A = 6 A Ad where A An is the Doppler 
width (including thermal and microturbulent contributions) 
at Te — 8 x 10 K. This form of truncated profile was used by 
Harper (1992) in the modelling of the 'hybrid' giant a TrA 
and variations on this approximate a pproach to PRD have 
been used by many authors (including McMurry 1999 in the 



modelling of a Ta u). The second approximate technique is 



that suggested by |Gayley (1998) 



which will be referred to 
as case (f). In that paper, the author suggests that if an ap- 
proximate PRD method that only requires CRD calculations 
is required, rather than simply truncating the Voigt profile, 
one should use a Voigt profile with a depth dependent Voigt 
parameter, a c ft = 1.8/r where r is the mean optical depth 
in the line. This approximation is intended to mimic the es- 
cape probability for first resonance lines, such as Lyman a, 
and is not strictly applicable to Lyman /3 (because of the 
influence of Balmer a), but it has been applied to both lines 
for comparison. The profiles calculated using these approx- 
imate PRD methods are compared with the case (d) PRD 
profiles for Lyma n a and Lyman fi in Figures ^ and q. 

As noted by Gayley (1998), neither of these approxi- 
mate methods is expected to make accurate predictions for 
the line profiles and this is evident in Figures ^ and |E| For all 
six calculated line profiles both the approximate techniques 
fail to accurately mimic the PRD profile. Owing to the trun- 
cation of the absorption profile, case (e) always significantly 
underpredicts the flux beyond the truncation frequency. For 
Lyman a, case (f) underpredicts the strength of the inner 
line wings and overestimates the strength of the outer line 
wings, producing a profile that is broader than the case (d) 
profile. 



4.2 Integrated Line Fluxes 

Table |^ shows the calculated integrated line fluxes at the 
stellar surface for each of the three models in all six cases 
considered. These fluxes show that not only the line profile 
but also the total flux is significantly affected by a PRD 
treatment; the total flux shows differences of up to a factor 
of 4.7 between CRD and PRD calculations. 

The main application of the approximate PRD tech- 
niques is to calculate total line fluxes (as discussed above 
it is known that they do not predict PRD line profiles cor- 
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Figure 4. Calculated Lyman a profiles. The solid line is case 
(d), the dotted line is case (e) and the dashed line is case (f). All 
fluxes are in ergs cm _2 s _1 A _1 at the stellar surface. AA is the 
wavelength relative to line centre in A. 



Figure 5. Calculated Lyman /3 profiles. The solid line is case 
(d), the dotted line is case (e) and the dashed line is case (f). All 
fluxes are in ergs cm _2 s _1 A _1 at the stellar surface. AA is the 
wavelength relative to line centre in A. 
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Table 3. Calculated integrated line fluxes from the models. All 
fluxes are at the stellar surface in ergs cm -2 s _1 . 



Line 


Case 


Cl 


Tau 


P Gem 


Procyon 


Lyman a 


(a) 


8.35 


X 


10 4 


7.26 x 


10 4 


9.67 x 10 5 




(b) 


1.78 


X 


10 4 


2.06 x 


10 4 


4.06 x 10 5 




(c) 


1.78 


X 


10 4 


2.06 x 


10 4 


4.06 x 10 s 




(d) 


1.77 


X 


10 4 


2.06 x 


10 4 


4.07 x 10 5 




(e) 


4.70 


X 


10 3 


1.69 x 


10 4 


3.86 x 10 5 




(f) 


5.67 


X 


10 s 


1.72 x 


If) 4 


4.19 x 10 5 


Lyman /3 


(a) 


2.15 


X 


10 1 


1.99 x 


10 2 


7.24 x 10 3 




(b) 


1.54 


X 


10 2 


2.31 x 


If) 2 


7.68 x 10 3 




(c) 


7.47 


X 


10 1 


2.00 x 


10 2 


7.12 x 10 3 




(d) 


7.19 


X 


10 1 


2.48 x 


10 2 


7.56 x 10 3 




(e) 


2.36 


X 


10 1 


2.06 x 


10 2 


6.86 x 10 3 




(f) 


9.67 


X 


10° 


1.92 x 


10 2 


6.65 x 10 3 



rectly). Comparing the PRD integrated line fluxes (case (d)) 
to those calculated using the approximate methods (cases 

(e) and (f)) shows that both approximate methods consis- 
tently underpredict the integrated line flux (the only excep- 
tion being Lyman a in Procyon case (f)). The approximate 
methods work best in Procyon (where the line fluxes agree 
with the PRD calculations to within 15 per cent) but the 
discrepancy is larger in the lower gravity stars: in a Tau, 
case (e) underpredicts the Lyman a flux by a factor of 3.8 
and the Lyman /3 flux by a factor of 3. It is not surpris- 
ing that these approximate methods, which are based on 
CRD, work best in Procyon where PRD effects are least 
important. The fluxes given in Table ^ indicate that for es- 
timating the Lyman a integrated line flux, case (f) is better 
than case (e), while for estimating the integrated Lyman j3 
flux, case (e) is better (this is not surprising since, as dis- 
cussed above, case (f) is not appropriate for Lyman 0). It 
is concluded that these approximate techniques can be use- 
ful in simulating PRD effects in the higher gravity stars: 
for Procyon, the case (f) Lyman a flux agrees with case (d) 
to within 3 per cent and for j3 Gem, the case (f) Lyman a 
flux agrees with case (d) to within 20 per cent, compared to 
a discrepancy between the flux predicted by normal CRD 
and case (d) PRD of a factor of 2 in Procyon and a fac- 
tor of 3.5 in p Gem. The approximate techniques appear to 
be less reliable in the lower gravity star a Tau where case 

(f) underpredicts the Lyman a flux by a factor of 3.12 and 
case (e) underpredicts the Lyman /3 flux by a factor of 3.05. 
However, as mentioned in Section 3, the significance of the 
integrated line fluxes in the a Tau model is limited by the 
reliability of the estimated line core profile produced by the 
extrapolation of the source function above the top point in 
the model atmosphere. 

Differences between CRD and PR D integrated line 
fluxes have been shown by ludge (1990)| to be related to the 
line thermalization depth. In that paper the author shows 
that in a solar model the Mg II fc-line (which forms in the 
chromosphere) thermalizes significantly higher in the chro- 
mosphere when using a PRD rather than CRD approach, 
while in a model of a Tau, the same line thermalizes, for 
both PRD and CRD calculations, not in the chromosphere 
but in the photosphere. Because the effective optical thick- 



ness of the line depends upon the thermalization depth, this 
means that the discrepancy between CRD and PRD calcula- 
tions of the Mg II k integrated line flux is significantly greater 
in the effectively thick solar model than the effectively thin 
a Tau model - the more effectively thin the atmosphere the 
less important are the details of radiative transfer in calcu- 
lating the emergent flux. Applying the same arguments to 
Lyman a would suggest that one should expect the discrep- 
ancy between computed CRD and PRD line fluxes to be 
significant only when the atmosphere is effectively thick in 
Lyman a (i.e. when the Lyman a line thermalizes high in the 
atmosphere) . Figure [] shows the computed source functions 
(S) and the Planck function (B) for each of the atmospheric 
models. The CRD source function is shown along with the 
monochromatic source function for case (d) PRD at four dif- 
ferent wavelength detunings (AA/AAd = 0,6,18,35 where 
AArj is the thermal Doppler width at 8 x 10 3 K). It can be 
seen that thermalization of Lyman a (in the sense that S 
does not significantly deviate from B) occurs in the chro- 
mosphere in each of the models and there is not a strong 
variation of thermalization depth with surface gravity. This 
difference between the thermalization properties of Lyman a 
and the Mg II line is due to the different photon destruction 
processes which thermalize the lines. In Mg II the line is ther- 
malized by electron collisional de-excitation and since this 
process occurs at a rate proportional to the electron density 
the photon destruction rate drops rapidly in the outer at- 
mosphere. This means that the source function can deviate 
significantly from the Planck function above a certain level 
in the atmosphere where the electron density has dropped 
low enough that collisional de-excitation can no longer ther- 
malize the line. Since the value of the electron density is 
lower in lower gravity stars this decoupling of the source 
and Planck functions will occur sooner (i.e. deeper in the 
atmosphere) in lower gravity stars. In Lyman a, collisional 
de-excitation is a less important thermalization process that 
photoionization of n — 2 by the Balmer continuum. Unlike 
the collisional de-excitation rate, the photoionization rate 
does not drop rapidly in the outer atmosphere and so this 
process is effective in preventing a substantial deviation of 
S from B until the high chromosphere where the thermody- 
namic properties of the atmosphere vary sufficiently rapidly 
that the photoionization process is no longer able to couple 
S to the local value of B. Thus, Lyman a is not effectively 
thin in any of the stars considered here and so radiative 
transfer effects are important in the calculation of both the 
line profile and the line flux in all the models. 

The integrated line flux for Lyman using CRD agrees 
with the PRD flux to within 5 per cent for Procyon and 25 
per cent for j3 Gem, significantly closer agreement than was 
found for Lyman a. In a Tau the difference between CRD 
and PRD fluxes is somewhat larger (factor of 3.3), but still 
less than the difference that was found in Lyman a. Fig- 
ure ^ shows the source functions and Planck function for 
Lyman 0. It can be seen that, like Lyman a, Lyman (3 ther- 
malizes in the high chromosphere in each of the models and 
so this line is not effectively thin in any of the models and 
radiative transfer effects will be important, not only in the 
calculation of the line profile but also in the calculation of 
the total flux. The fact that, despite the Lyman j3 line be- 
ing effectively thick, PRD and CRD calculations predict line 
fluxes that agree substantially better than those found for 
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Lyman a is a strong indication of how much less important 
PRD effects are in Lyman than Lyman a, however it must 
be stressed that the effects of PRD on the Lyman f3 flux are 
still significant, particularly in the low gravity stars. 

4.3 Implications for Modelling 

The calculations discussed above were all carried out us- 
ing fixed atmospheric models. However, this process is not 
strictly self-consistent. It is known that using a PRD (rather 
than CRD) tr eatment alters the level populations of the try 



drogen a tom ( rlubeny fc Lites 1995 , Luttermoser fc John-, 
son 1992j ). Figure § shows the level populations calculated 



for the four different redistribution cases using the /3 Gem 
atmospheric model. The variation in the level pop ulations is 



the n 



similar to that found by Hubeny fc Lites (1995) : 
and n = 3 level populations rise because of the increased 
trapping of photons in the Lyman a line (as described in 
Section 4.1). An increased n — 2 population means that the 
proton population also rises because the primary ionization 
mechanism is photoionization from n — 2. These increases 
in the n — 2, n = 3 and proton populations are balances by 
a decrease in the n = 1 ground state population. Physically, 
a significant increase in the proton density carries with it an 
increase in the electron density and so to hold the electron 
density constant (at a value calculated from the hydrostatic 
equilibrium assumption in CRD) is not acceptable. There- 
fore, it is necessary to use a PRD treatment of hydrogen 
when constructing a model atmosphere in order that PRD 
calculation results can be consistent with the assumption of 
hydrostatic equilibrium. 

For this reason, each of the atmospheric models has 
been recomputed using a PRD treatment of hydrogen in the 
hydrostatic equilibrium integrations. For a Tau, the original 
model has simply been recomputed in hydrostatic equilib- 
rium with PRD. This means that the new model of a Tau is 
no longer strictly consistent with the emission measure dis- 
tribution from which it was derived by McMurry (1999) The 
models of /3 Gem and Procyon have been recomputed in hy- 
drostatic equilibrium with PRD and have also been iterated 
to consistency with the emission measure distributions from 
which they were derived. The only difference between the 
three new PRD-based models and the original three models 
is in the electron density. Figure |^ shows the change in the 
electron density in each model. f3 Gem shows the largest dif- 
ferences, up to a factor of three between the electron density 
in CRD and PRD calculations. As would be expected, these 
differences occur in the chromosphere and low transition re- 
gion, where the Lyman lines are forming. Procyon shows 
the smallest differences in the electron density of not much 
more that 20 per cent, a Tau shows significant differences 
between the CRD and PRD calculations, in a similar way to 
P Gem, but the maximum change in a Tau is only about a 
factor of two. It may seem surprising that the lower gravity 
star a Tau shows smaller changes than /3 Gem, but it must 
be borne in mind that the direct comparison of these two 
models is not valid: the original a Tau model is not a true 
CRD model but was constructed by McMurry (1999) using 
an approximate PRD method. 

Figure ^ makes clear that a PRD treatment of hydro- 
gen is important in constructing atmospheric models, par- 
ticularly but not exclusively, for low gravity stars. In, for 
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Figure 6. Source functions for Lyman a. In each case the dotted 
line is for CRD, the solid line is for PRD (AA = 0), the dashed line 
is for PRD (AA = 6AA D ), the dash-dot line is for PRD (AA = 
18AA D ), the dash-triple-dot line is for PRD (AA = 35AA D ) and 
the long-dash line is the Planck function. All source functions 
are in ergs cm -2 s — 1 sr _1 Hz -1 and column mass density is in 
g cm -2 . AA D is the thermal Doppler width at T e = 8 X 10 3 K. 
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Figure 7. Source functions for Lyman (i. In each case the dotted 
line is for CRD, the solid line is for PRD (AA = 0), the dashed line 
is for PRD (AA = 6AA D ), the dash-dot line is for PRD (AA = 
18AA D ), the dash-triple-dot line is for PRD (AA = 35AA C ) and 
the long-dash line is the Planck function. All source functions 
are in ergs cm~ 2 s _1 sr _1 Hz -1 and column mass density is in 
g cm~ 2 . AA D is the thermal Doppler width at T e = 8 X 10 3 K. 
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Figure 8. Calculated hydrogen level populations for (i Gem. The 
solid line is case (d) PRD and the dotted line is CRD. Popula- 
tion number densities are in cm~ 3 and column mass density is in 
g cm -2 . 
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Table 4. Calculated integrated line fluxes from the four new 
PRD-based models with case (d) PRD treatment of hydrogen. 
All fluxes are at stellar surface in ergs cm -2 s _1 . 
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Figure 9. Electron densities from hydrostatic equilibrium cal- 
culations. The solid lines are from the original (CRD) models 
and the broken lines are from models constructed based on hy- 
drostatic equilibrium models computed using a PRD treatment 
of hydrogen. The electron number density is in cm" 3 and the 
column mass density is in g cm -2 . 



Transition 


a Tau 


(3 Gem 


Procyon 


Lyman a 


1.23 X 10 4 


3.20 x 10 4 


4.03 X 10 5 


Lyman (3 


3.79 x 10 1 


3.44 x 10 2 


7.50 x 10 3 



example, the model of j3 Gem, the electron density in the 
chromosphere is very significantly different from that pre- 
dicted by the CRD treatment and it will lead to substantially 
different fluxes of collisionally excited lines of all elements. 
It is not within the scope of this paper to give a detailed dis- 
cussion of the atmospheric emission lines of elements other 
than hydrogen, but the importance of the PRD treatment 
of hydrogen on these lines will be discussed in subsequent 
papers. In particular the collisionally formed lines of C II 
and the fluorescent lines of O I (pumped by Lyman /3) will 
be of interest. Table ^ gives the integrated fluxes of the Ly- 
man a and Lyman /3 lines calculated using the full case (d) 
PRD with the new PRD-based models. If these fluxes are 
compared with the case (d) results in Table | it can be seen 
that the difference in the fluxes reflects the difference in the 
electron density: f3 Gem shows a very significant difference, 
a Tau shows a significant difference but Procyon shows very 
little difference. 



4.4 Comparison with Observations 

The work discussed in the previous sections has led to the 
conclusion that, for a given atmospheric model, predictions 
based on a CRD treatment of hydrogen differ from those 
based on a PRD treatment in three ways. First and foremost, 
CRD leads to predicted line profiles for Lyman a that have 
substantially stronger wings (and hence greater integrated 
line flux) than the more satisfactory PRD treatment. PRD 
also predicts significantly different profiles for Lyman /3 in 
low gravity stars - this effect should be observable in future 
observations of late-type stars with the FUSE instrument. 
Secondly, because the line profiles are different, there may 
be differences in the properties of lines that are excited by 
fluorescence in the Lyman lines, depending on the wave- 
length at which lines are pumped. Thirdly, because a PRD 
treatment leads to a different chromospheric electron den- 
sity there will be a change in the strength of collisionally ex- 
cited chromospheric/transition region lines (such as the C II 
lines at around 1335 A). The secondary effects on the lines 
of other atomic species will not be discussed here but will 
be addressed in subsequent papers on the details of the at- 
mospheric models. The most easily observed consequence of 
PRD (as opposed to CRD) is the significant relative reduc- 
tion of the line wing intensity in Lyman a compared to the 
line core. Figure ^ show observed Lyman a profiles and the 
synthetic profiles calculated using the original models with 
case (a) CRD, and also the new PRD models treated with 
full case (d) PRD. The observed profile for a Tau is from the 
IUE high dispersion exposure SWP6679 (C. Jordan (PI)). 
The observations for Procyon and f3 Gem are from GHRS 
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Table 5. Parameters used to correct calculated line profiles for 
interstellar H I absorption and macroturbulence. For (3 Gem and 
Procyon there are two separate components to the interstellar 
absorption. 



Star 


Parameter 


Value 


a Tau 


Cmacro (km S" 1 ) 


17 




f micro (km s" 1 ) 


13 




N H (cm- 2 ) 


5 x 10 18a 




6 (km s- 1 ) 


ll b 




uo(km s — 1 ) 


- 30 c 


P Gem 


£macro (km S ) 


16 




£micro (km S" 1 ) 


12 




N H (cm" 2 ) 


1.15 x 10 18 /6.82 x 10 17d 




b (km s- 1 ) 


12. 33/11. 02 d 




vq (km s _1 ) 


21.9/33.0 d 


Procyon 


£macro (km S -1 ) 


29 




£micro (km S" 1 ) 


14 




N H (cm" 2 ) 


7.5 x 10 17 /4.0 x 10 17e 




b (km s" 1 ) 


10.78/10.78 e 




vo (km s - *) 


20.5/21. e 



a McClintock et al. (197 8) 
6 assumed 



~irro.A on Mg TT k- 



Dring ct al. (19971 



ne (see text) 



Linsky ct al. (1995) 



exposures Z17X0303M (J. Linsky (PI)) and Z2SI0406T (R. 
Henry (PI)) respectively. 

The atmospheric models contain a microturbulent ve- 
locity (£ m icro) which is taken into account in the radiative 
transfer calculations. In order to correctly match the ob- 
served widths of optically thin lines, computed profiles gen- 
erally have to be broadened by convolution with a Gaussian 
profile whose width corresponds to a most probably macro- 
turbulent speed £ m acro. Macroturbulence is interpreted as 
large scale turbulence in the stellar atmosphere. The values 
of £macro which have been used are shown (along with the 
microturbulent velocity £ m i C ro at T e = 10 4 K) in Table || 
These values of ^ macrn have been deduced from the widths 
of optically thin emission lines that form in similar regions 
of the atmosphere to the Lyman a line and will be discussed 
further in later work regarding the atmospheric models. The 
computed profiles have been corrected for interstellar ab- 
sorption. Th is has been pe rformed using the optical depth 
described by Munch (1968) in terms of the hydrogen column 



density Nh , the line width parameter b (which in the current 
work is taken as V2 times the Doppler dispersion parame- 
ter) and the wavelength shift relative to line centre vo- The 
adopted parameters for the interstellar absorption are shown 
in T able ^. For Procyon and [3 Gem, Linsky et al. (1995) 
and Dring et al. (1997) respectively have deduced param- 
eters for the interstellar absorption using two component 
models. Both component s of the H I absorption h ave been 
included here. For a Tau McClintock et al. (1978) give only 



the column density and so the line width has been taken as 
a typical value of 11 km s _1 . The velocity of the interstellar 
component has been taken to be the same as the velocity of 
the interstellar absorption component reported in the a Tau 
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Figure 10. Observed and calculated Lyman a profiles. The solid 
line is the observed flux, the dotted line is the calculated flux 
using CRD and the dashed line is the calculated flux using PRD, 
in units of 10 -11 ergs cm" 2 s~ 1 A -1 . The wavelength is in A. The 
central emission feature in a Tau and p Gem is due to geocoronal 
emission. 
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Mg II fc-line by Robinson, Carpenter fc Brown (1998)| . For 
comparison with the observations, the final a Tau profile 
has been convolved with a Gaussian of FWHM= 25 km s _1 
to take the resolution of the IUE instrument into account. 

It is clear from Figure |Io| that the PRD line profiles are 
in much better agreement with the observed line shapes than 
are CRD profiles, a clear indication that a PRD treatment is 
very important to successfully modelling the Lyman a line. 
The calculated line profiles do not agree closely with the 
observed profiles, but this is most likely due to inadequacies 
in the model atmospheres. In particular the total Lyman a 
flux is being significantly underpredicted in /3 Gem. It is 
beyond the scope of this paper to recompute the atmospheric 
models, but this will be done in future work. The asymmetry 
of the observed Procyon and (3 Gem profiles is well accounted 
for by the effects of interstellar absorption. The significant 
asymmetry of the a Tau profile is likely to be due to flows 
in the atmosphere (or stellar wind) which are not included 
in the hydrostatic model used here. The observed profiles 
also show additional features due to deuterium absorption 
(which is seen in the blue wing of the Procyon and j3 Gem 
profiles at AA = — 0.3A from line center) and geocoronal 
emission (which is seen as a tall spike near the centre of the 
f3 Gem and a Tau profiles). 



5 CONCLUSIONS 

It has been demonstrated that there are substantial differ- 
ences between the results of PRD and CRD calculations of 
the Lyman lines, and that these differences are significant 
in the modelling of stellar atmospheres. The importance of 
PRD effects grows in lower surface gravity stars, but the 
effects are not negligible in any of the evolved stars being 
studied here. The primary effect has been shown to be a 
significant reduction in the Lyman a wing intensity which 
cause a major reduction in the computed Lyman a flux and 
has significant repercussions on the atmospheric structure 
because of the way in which it modifies the ionization bal- 
ance in hydrogen. The differences in Lyman (3 are less sig- 
nificant and it appears that for lower gravity stars (such as 
Procyon) this line may be modelled reasonably well with a 
CRD assumption. However, it has also been demonstrated 
that this is not always the case: the differences in a Tau are 
substantial. 

The significance of these results goes beyond simply pre- 
dicting the Lyman line profiles and fluxes. Because Lyman a 
is such an important radiative energy loss mechanism, its 
behaviour has repercussions for the understanding of chro- 
mospheric heating. The treatment of this line has the power 
to change the structure deduced for the atmosphere, caus- 
ing a significant change in the electron density which will 
manifest itself in the predicted fluxes of collisionally excited 
emission lines. The issue of fluorescence is also significant 
and, depending on the wavelength of the pumped lines, the 
differences between the PRD and CRD profiles could lead 
to different predicted fluxes of the fluorescent lines. This is 
important since the lowest gravity stars show the richest 
spectra of fluorescent lines and it is these stars which show 
the strongest PRD effects. 

The atmospheric models that have been employed as 
test cases in this paper are preliminary (it is clear from 



the comparison in the previous section that they do not, at 
present, entirely reproduce the observations) but these mod- 
els provide a useful set of realistic test cases to examine the 
importance of PRD in hydrogen. The further development 
and discussion of the atmospheric models will be presented 
in future work. 
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